www.gusucode.com > 网格法与lingo,内含例子及多个matlab代码 > 网格法与lingo,内含例子及多个matlab代码/third day/f2.m

    function y=f2(x1,x2,x3,x4,x5,x6,x7)
y=-148257/1000*x1/x5/(x3/(x2-x1))^(3/20)*((1-131/50*(1-9/25/(x4/x2)^(14/25))^(3/2)*(x4/x2)^(29/25))/x6/x7)^(1/2)*x3/(x2-x1)^2+8721/100*x1/x5*(x3/(x2-x1))^(17/20)/((1-131/50*(1-9/25/(x4/x2)^(14/25))^(3/2)*(x4/x2)^(29/25))/x6/x7)^(1/2)*(24759/31250*(1-9/25/(x4/x2)^(14/25))^(1/2)/(x4/x2)^(2/5)*x4/x2^2+3799/1250*(1-9/25/(x4/x2)^(14/25))^(3/2)*(x4/x2)^(4/25)*x4/x2^2)/x6/x7;